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Abstract 

We consider a simple model to describe the widths of the mode locked inter- 
vals for the critical circle map. Using two different partitions of the rational 
numbers, based on Farey series and Farey tree levels respectively, we calculate 
the free energy analytically at selected points for each partition. It is found 
that the result of the calculation depends on the method of partition. An 
implication of this is that the generalized dimensions D q are different for each 
partition except when q = 0, i.e. only the Hausdorff dimension is the same in 
each case. 

PACS: 37-XX,37E10,37D35,37L30 
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1. Introduction 

On the basis of their numerical investigations, Jensen et al. [9, 10] con- 
jectured that for the critical circle map the mode locked intervals form a 
complete Devil's staircase whose complementary set is a Cantor set with a 
fractal dimension D ~ 0.87. Further, their numerical results indicated that 
this dimension is universal for maps with a cubic inflection point. 

Their results were obtained in a purely empirical fashion by computing 
the width of smaller and smaller mode locked steps and subtracting the 
summed widths from the total staircase length. 
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The set in question is the set of irrational windings, which do not lead 
to mode locking. It is multifractal in nature and, more precisely, is believed 
to have a universal Hausdorff dimension Djj ~ 0.87, In order to better 
characterize such sets, different thermodynamic formalisms [6, 5] have been 
introduced and with them, the idea of generalized dimensions D q [8]. 

In this paper we point out that for the circle map, such a description 
is not unique and depends upon the details of how the rational numbers 
are partitioned. Although we find the thermodynamics of the set for two 
distinct partitions is quite different, the Hausdorff dimension is the same in 
both cases. 

2. Model for Interval Widths 

In order to illustrate this ambiquity, we shall examine a simple model for 
the widths of the mode locked intervals of the circle map and show that two 
different partitions of the rationals lead to two different functional forms for 
the free energy F(f3) (defined below) and hence to different D q . Suppose 
that the width A(P/Q) of the mode locked interval for the critical circle 
map associated with the winding number P/Q is given by [11] 

A(P/Q) = kQ- 5 , (1) 

independent of P. P and Q are coprime integers; k and 5 are constants (we 
assume 5 > 2 ). 

In this model we know the widths of the mode locked intervals, however 
we do not know precisely where the left and right hand ends of the interval 
are located on [0,1]. This means we do not know exactly where the holes 
are located. However, when experiments are performed on one-dimensional 
quasiperiodic systems at the borderline of chaos [2], it is only possible to 
determine the 'devil's staircase' for mode locked intervals at rational winding 
numbers P/Q where Q is relatively small. Nevertheless, it is possible to 
determine the thermodynamics of the system including an estimate of the 
Hausdorff dimension of the set of irrational windings. 

In this paper, we use a simple model for these widths and are able to 
construct the thermodynamics of this model. This leads to an estimate of 
the Hausdorff dimension in agreement with that obtained by more general 
considerations. 

From Eq (1) we may analytically calculate the total sum of the widths 
of the mode locked intervals (which is taken to be one at criticality) . 

YJl A ( p /Q) = k E <t>(Q)Q~ 5 = k ^7sT = x ' (2) 

Q P Q=l ^ ' 

where ( is the Riemann zeta function and <p(Q) is the Euler function [7]. The 
sum in Eq (2) is over all irreducible rationals P/Q, P < Q, and A(P/Q) is 
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the width of the parameter interval for which the iterates of a critical circle 
map lock onto a cycle of length Q, with winding number P/Q. Therefore k 
is fixed to be 

Now the Hausdorff dimension is the largest value of [3 for which the sum 

Q P ^ 



= kPJ2<t>(Q)Q- sp (4) 



C(5P-l) 



cm 

diverges. It is clear from Eq (4) that there is a pole singularity in the sum 
when ((5/3 — 1) — > 1 from above, i.e. (3 — > 2/5. We conclude therefore that 
the Hausdorff dimension is given by 

0h = 2/5. (5) 



3. Thermodynamic Formalism 



A method for extracting a spectrum of scalings from experimentally or 
numerically generated strange (fractal) sets was introduced by Halsey et 
al. [6] using a "thermodynamic formalism". They found that the scaling 
properties of normalized distributions lying upon such sets were character- 
ized by two indices, namely a which determines the strength of the singular- 
ities and / which describes how densely they are distributed. The spectrum 
of singularities is described by giving the range of a-values and the function 
/ = /(<*). 

Shortly thereafter Feigenbaum [5] pointed out that the work of Halsey et 
al. was in fact the microcanonical version of a thermodynamical formalism 
introduced by Vul, Sinai and Khanin [13] in 1984 some years earlier. See 
also Ruelle [12]. He elucidated the canonical version (canonical paradigm = 
CP) in contrast to the microcanonical version (microcanonical paradigm = 
MP). In this formulation, the probability measure is assumed constant on 
the set. This formalism introduces a free energy F{(5) which is related to the 
f(a) as indicated below. In the following we use the notation of Feigenbaum 
[5]. 

Consider a dynamical system whose attractor can be hierarchically rep- 
resented as a set of N n intervals l\ n ^ , i = 1 . . . N n of lengths A at the n-th 
level. The free energy is defined by 

AT n -^) = £|AfV (6) 
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where asymptotically F becomes independent of n. The relation between 
the MP quantities ( Halsey et al. [6]) and the CP quantities (Feigenbaum 
[5]) is given by 



a = 1/F'((3) 
f = (3-F((3)/F'((3) 
q = -F((3) 
t = -0 
D q = (3/[l + F((3)] 



(7a) 
(7b) 
(7c) 
(7d) 
(7e) 



3.1. Free Energy for Farey Series Partition In order to proceed 
further, we first recall that for a given Q, the number of irreducible rationals 
P/Q is 4>(Q). The Farey series of order Q, {Fq}, is a set containing the 
monotonically increasing sequence of all irreducible rationals P'/Q' (P' and 
Q' are coprime) between and 1 whose denominator does not exceed Q [7]. 
Next we define a sequence of integers {Q n } such that asymptotically 



With these Q n , we define a sequence of Farey series {FQ nl }, {-Pq„}, {-^Q n +i} 
etc. It follows that the number of rationals contained in {Fq u+2 } but not in 
{FQ n+1 }, is asymptotically twice the number contained in {FQ n+1 } but not 
in {Fq u }. We therefore have a simple prescription for grouping the rationals 
into distinct exclusive sets {{FQ n+1 } - {Fq u }} whose membership increases 
exponentially fast. We may re-express this statement by saying that the 
ratio of the number of rationals contained in each of these consecutive sets 
is given by 



Qn+l 
Qn 



(8) 




Qn + 2 



Q 



Q 



Qn 



Qn+2 Qn+l 
Qn+1 ~ Qn 
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asymptotically. We have used [7] 



J>(Q)~4;Q : 



Q 



in the large Q limit. 
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This construction defines a partition of the rationals from which one may 
define the free energy F{(5) in this simple model asymptotically through 



2~ F (l 3 ) = ® n+1 

Qn+l 



E <P(Q)Q~ sp 

3n+l 
2n + l 

E <KW _a/3 



rWn + 2 
'Qn+l 



L 



•Qn+l 
Qn 

n 2-8f3 _ n 2-8f3 

~~ n 2-<5/3 _ n 2-8f3 
^n+l 



from which we conclude that 



m = f - 1- (9) 

The free energy has the same form as that for a single scale Cantor set, 
which is not surprising since there only appears to be a single length scale 
associated with this model. The Hausdorff dimension is found from the 
condition 

F{0 = Ph) = 

i.e. 

in agreement with the result of Eq (7d). Using the f(a) language 

2 

a = /(") = g 

i.e. the /(a) curve collapses to a single point. It should be noted that the 
expression for F(/3) quoted in Eq (9) in no way depends on the partition 
defined through Eq (8). Indeed we could choose any a > 1 such that 

Qn+l r- 

— — ~ Va, 
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which would lead instead to 



a 



and again we have 



A natural choice for the Q n that suggests itself is the Fibonacci sequence, 
Qn = F n in which case 

Qn+l 1 



where p is the golden mean 



and 



Qn P 
y/E- 1 



1 

~3 ' 



In fact such a choice has been used [1] for evaluating the free energy (or its 
equivalent) and the Hausdorff dimension in the case of the real mode locking 
critical circle map. This is a different partition of the rationals which also 
utilized Farey series. The crucial point, however, in both examples is that 
one may generate partitions of the rationals whose membership increases 
exponentially fast from one partition to the next; by construction the vari- 
ation in size of the denominator Q in the set {{FQ n+1 } — {FQ n }} is only of 
order 1. By contrast the Farey level construction which we shall discuss next 
allows enormous variations in Q at a given level. The smallest denominator 
is Q while the largest denominator is of order (1/ p)Q . 

3.2. Free Energy for Farey Tree Partition The Farey tree level 
partition of the rationals has been frequently used in discussions of the circle 
map [1, 3] and one may readily calculate the free energy using this scheme. 
Since the number of rationals doubles from one level to the next, the free 
energy may be defined through 



Era 



-5(3 



2"™ = '—. Ml - M (10) 



when comparing level n of the Farey tree with level (n + 1). The (resp. 

Q^ n+1 ^) are the denominators of the rationals pj n ^ /Q^ (resp. P- n+1 ^ / Q[ n+1 ^ ) 
at level n (resp. (n + 1)). In general, it is not possible to evaluate Eq (10) 
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analytically and solve for F{(3). However it is possible to evaluate it exactly 
at certain values of 5(3. From Eq (10) we have 

2-FWS) = _i (n) 



4. Comparison of Free Energies for Different Partitions 

Eq (11) may be exactly evaluated for a number of values of (3. At [3 = 
we have 

on+l 

= 2 

i.e. 

F(0) = -1. 

Here we have used the fact that the number of rationals at level n in the 
Farey tree is 2 n . This result for F(0) is in agreement with Eq (9) when (3 is 
set = 0. 

When (3 = — 1 we have 

13 = -I: 2' F ^ 



3 

3 n+l 



3 n 

= 3 



from which we conclude that 



Here we have used the fact that [1] 

j 

at the level n of the Farey tree. However, Eq (9) predicts that 

in disagreement with Eq (13). Therefore we see that we get different values 
for the free energy at (3 = — 1 according to the partition. 
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It is also possible to evaluate F((3) exactly at another point in the Farey 
tree construction 

E lQl n+1) f 

P = -3: 2- F HV 5 ) =-! — - T - 

E Wf 

so 

F(-h = - l ^l. (14) 
1 5' log2 1 ; 

Here we have used the fact that [1] 

£ [Q/ B >] 3 ~ 7". 
i 

The original partition predicts from Eq (9) that 

again in disagreement with Eq (14). These results are summarized in Table 
1. 



Table 1 near here 



Therefore it is possible to show exactly analytically that a different par- 
tition of the rationals yields different values for the free energy at specific 
values of (3 (and hence a different singularity spectrum /(en)) in this simple 
model for the widths of mode locked intervals. From Eq (7e) we may express 
this in terms of the generalized dimensions D q . For the first partition 

D q = P/(1 + F(P)) 
2 

6 

independent of q. Equivalently in terms of MP, Eq (9) states 

r-§ta-l). 

Since [6] 

T=(q- 1) D q , 

we have 




once more. 
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We see that D q is no longer constant. Cvitanovic [4] has found the 
q(r)(= —F((3)) curve for the case 5 = 2 using the Farey tree construction 
together with accelerated convergence numerical techniques. It is clear that 
in this case q(r) is no longer linear in r suggesting that there is more than 
one length scale which at first sight seems surprising in this simple model. 
Since q(r) is no longer linear in r (as it was in the Farey series partition), 
D q is no longer a constant but varies with q. 



Figure 1 near here 



In Figure 1, we plot F((3) over a range of (3 values for both partitions 
of the rationals. The straight line F((3) = ((3 — 1) (shown as broken) is 
associated with the Farey series partition. The other line (continuous) shows 
F((3) computed using Farey tree levels. For the purpose of illustration, the 
value 5 = 2 has been chosen. For negative half integer (3, F(f3) is known 
analytically for the latter partition [1]. See Table 2. For positive (3, the 
thermodynamic function is computed for Farey tree levels 13, 15, 17. Using 
accelerated convergence techniques one then obtains F((3) = for (3 > 1 
[1, 4]. Clearly the thermodynamic functions are quite different and this 
conclusively demonstrates that their determination depends strongly on the 
method used to partition the rationals. 



Table 2 near here 



5. Conclusion 

In conclusion we see that by analyzing an oversimplified model for the 
widths of the mode locked intervals in a critical circle map and partitioning 
the rationals P/Q in two distinct ways, it is possible to derive two completely 
different functional forms for F(0) [resp. q{r)] which only agree at two 
different values of (3 [resp. r]. As noted one of these is (3 = (which was 
shown analytically in Eq (12)) while the other is presumably at — r = (3 = 
(3h- Cvitanovic [4] has shown that the latter is true for 5 = 2 to a high 
precision using numerical techniques. Even if one has doubts about such 
numerical techniques, the ability to compute the thermodynamic function 
analytically at isolated points and find that the result depends on the form 
of the partition used demonstrates unequivocally that a serious ambiguity 
has arisen. One would expect that other partitions different from those 
discussed here would yield thermodynamic functions that are different from 
the ones depicted in Figure 1. 

An alternative way of stating the conclusion is that the Farey series 
scheme for partitioning the rationals when one adopts this simplified model 
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for the mode locking widths yields a single constant generalized dimension 
D q = Do, while the Farey level tree partition clearly yields a result for D q 
which is not constant with variable q. Thus there is clearly ambiguity in 
this attempt to determine the generalized dimensions associated with the 
circle map which presumably persists for the real critical circle map which 
is genuinely multifractal. 
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Figure 1. A comparison of the two partitions used to determine F((3). The continuous 
line is F((3) = (/3 — 1) and is associated with the Farey Series/ Sequence partition. The 
broken line shows F{f3) for the Farey tree partition. This maybe be calculated analytically 
for negative half integer /3 (see Table 1) and computationally when beta is positive for 
finite order Farey tree levels. See [1] who conjecture that asymptotically, F({3) converges 
to /3 = 1 so that F(P) for the Farey tree partition is C 1 at (3 = 1. These different partitions 
lead to different thermodynamic free energy F(f3) which have only two points in common. 
One is at F(0) = —1 where all partitions must agree. The other is at F(l) — which 
means both partitions lead to the same Hausdorff dimension Dh = Dq = 1. Note that for 
the purposes of illustration, the value 8 = 2 has been chosen. 
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p 


Farey series 


Farey tree 


1 











-1 


-1 


-1/2 


-3/2 


-log 3/ log 2 


-3/2 


-5/2 


- log 7/ log 2 



Table 1. Thermodynamic free energy F(f3) evaluated at selected values of f3 for both 
partitions with 8 = 2. 
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Asymptotic limit of 2 ^ ^ 


1 


1 





2 


-1/2 


3 


-1 


(5 + v / 17)/2 


-3/2 


7 


-2 


(ll + VTl3)/2 


-5/2 


7 + 4^/6 


-3 


26.20249 . . . 


-7/2 


41.0183... 



Table 2. Analytic limits of the ratio ^ [Q ( j n+1) ] / ^ [Q< n) ] for consecutive levels 
in the Farey tree in the asymptotic limit n — > oo for various /3 (see Eq (10)). 
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